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Abstract: - The purpose of the paper is to substantiate a new approach to forecasting hydrocarbon potential 
of subsoil plots via co-analysis of satellite, geo-geophysical and geochemical information by means of 
subjective logic. Basic concepts of subjective logic are given, the method of forecasting the oil and gas 
potential of a subsoil plots is described. The method is tested by applying it to real area with hydrocarbon 
deposits. 
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1 Introduction The main reasons for the erroneous conclusions 
For over a century, hydrocarbons have remained are the insufficient knowledge of the genesis of 
one of the most common sources of energy, and hydrocarbon formation, the ambiguity of the 
every year humanity's demand for oil and natural existing forecast criteria for hydrocarbon 
gas 1s constantly growing. Despite increased use of accumulations and their dependence on geological 
renewable energy sources and increased energy conditions and occurrence depths, the low 
efficiency, demand for oil and gas continues to efficiency of existing direct methods for 
grow along with energy consumption in most hydrocarbon deposits prospecting, the presence of 
regions of the world. For example, the annual artifacts, etc. [3]. 
report of the International Association of Oil and Satellite technologies occupy an important place 
Gas Producers [1] states that the demand for oil in in hydrocarbon prospecting, since they give a 
2018 increased by 30% compared to 2000, and the possibility promptly to receive important 
demand for natural gas increased even more - by information about landcover and terrain in whole 
60% compared to the same period. [4, 5]. Thus, a procedure for predicting the oil and 
It should be said that the 2020 economic crisis gas potential of subsoil plots requires the utilization 
caused by the SARS-CoV-2 coronavirus has of a large amount of heterogeneous data that have 
greatly changed the situation in the world oil and different accuracy and completeness, differ in 
gas market, but it is obviously that this is a transient presentation forms, spatial scale, the moment of 
phenomenon. Understanding this fact and looking assessment or observation, etc. 
ahead, oil and gas producing countries do not Analysis of heterogeneous and incomplete data 
reduce, in general, the volumes of work on is a rather complicated process; therefore, to solve 
exploring new hydrocarbon deposits. the problem, artificial intelligence methods, 
However, the world experience shows, the algorithmic decision-making procedures, etc. are 
effectiveness of search of new hydrocarbon used [6]. 
deposits is still not high enough, and the positive The purpose of the paper is to develop a new 
conclusions according the availability of approach to forecasting the oil and gas potential of 
hydrocarbons which are provided by experts on the subsoil plots. This approach actively exploits the 
basis of studying the subsoil plots properties are mathematical apparatus of subjective logic 
unfortunately often not confirmed by the results of developed by A. Josang [7], and its implementation 
drilling [2]. allows you to: 
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- adjust the subjectivity and uncertainty of 
expert judgments, 

- generate adequate on the results of the analysis 
one-time or rare events; 

- avoid the majority of logical-computational 
problems that usually arise when combining 
conflicting hypothesis and opinions. 

In total, these advantages provide a good 
opportunity to increase the reliability of forecasting 
gas and oil potential of subsoil plots. 

The paper is structured as follows. Section 2 
examines state of the problem. Section 3 describes 
basic concepts of the subjective logic [7]. In section 
4, the new method to forecast gas and oil potential 
of subsoil plots through co-analysis of satellite, 
geological, geophysical and geochemical 
information by means of subjective logic is 
proposed and described in detail. In section 5, we 
conduct validation of the proposed method on real 
data. And, finally, in section 6, conclusions and 
directions for further research are represented. 


2 State of the problem 

The intellectualization of data processing is based 
on the creation and the implementation of 
mathematical (formal) models of objects and causal 
relationships between them. The formal models 
allow to algorithmize the hydrocarbons prospecting 
process and, as a consequence, to reduce, at least in 
part, the level of subjectivity of forecast estimates 
and make them more adequate and reliable. A 
confirmation of this fact is the successful 
application of expert systems PROSPECTOR, 
Dipmeter Advisor, GEOPLAY, but not only [8]. 

It is known that most of the known models and 
algorithms for data processing and decision- 
making are constructed using a probabilistic- 
statistical approach. However, in the tasks of 
remote forecasting of the oil and gas potential of 
subsoil plots, relying on a probabilistic paradigm 1s 
not the good solution. 

Firstly, in most cases of studying subsurface 
areas, the data obtained are at least partially 
inaccurate and ambiguous, and the amount of 
available information is limited. In such cases, it is 
difficult to talk about statistically reliable solutions. 

Secondly, despite significant progress in 
geological science, many issues of oil and gas 
genesis and of the influence of the deposit presence 
on soil and vegetation characteristics remain open 
and debatable. 

Due to these reasons, expert judgments about 
the oil and gas potential of the subsoil usually carry 
a significant element of uncertainty and are often 
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even contradictory. However, is known the 
processing of conflicting data by methods of 
probability theory is complex and can lead to 
incorrect results [9]. 

Therefore, new approaches’ to search, 
forecasting oil and gas potential of subsoil plots, 
based on other theoretical models are being 
developed. In the monography [10], a methodology 
for forecasting hydrocarbon potential of subsoil 
plots, based on the algorithmization of the 
processes of analysis of satellite, geological, 
geophysical and geochemical information is 
described. In [11], the method using the apparatus 
of fuzzy sets and fuzzy logic for of oil and gas 
prospecting sections of the Dnieper-Donetsk 
depression is developed. In [12], the technology of 
use of remotely-sensed data for quantitative 
estimate the potential of petroleum, natural gas and 
minerals taking into account the specific features of 
area is developed. In [13], it is proposed to detect 
intersection nodes of fault lines that are promising 
for hydrocarbon accumulation using the 
mathematical tools of the Dempster—Schafer 
evidence theory [14]. 

Authors propose a new approach for forecasting 
gas and oil potential of subsoil plots via co-analysis 
of satellite, geological, geophysical and 
geochemical information. The proposed approach 
relies on subjective logic, mathematic apparatus of 
which allows operating with uncertain and 
incomplete data, simple hypotheses and their 
disjunctions, as well as to combine conflicting 
hypotheses and opinions. 


3 Basic concepts 

Suppose @ is a set of mutually exclusive and 

collectively exhaustive hypotheses: 

© ={M, Hh, ... , Hj, ... Hy} (1) 
and denote as 2! a set, containing all subsets of 
different hypotheses: 

Or; = (5, {Hi}, {H2}, ea {Hj}, wee {Hy}, 

{Hı U Ap}, a {Hı AITE Dad Hy}, O), (2) 
where © denotes an empty set. 

The set © is called Frame of Discernment 
(FoD); elements H; | j = 1, 2,..., J of this set are 
called simple hypotheses. The set 2' is called 
Power Set; its subsets in the form of a disjunction 
of two or more hypotheses are called compound 


hypotheses. 

Each subset xIx e 2° is associated with a 
corresponding number mo(x), that satisfies 
conditions: 
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Me (x) 20 
Mm, (DO) =0 (3) 
S m(x)=1 
xe2'0 
The number mo(x) is called belief mass for the 
subset x. If moe(x) >O0 the subset x is called focal 
one. The set of belief masses for all elements of the 
set 2'°' forms the Belief Mass Assignment (BMA), 
which is indicated as me(Q). 
In subjective logic, for given xIx e 2" the 
following functions are introduced: 


belief function 
b(x) =$ mO); x,y e 2 (4) 
yox 
disbelief function 
d(x) = $ m (y); x,y € 2° (5) 
yAax=Ø 
uncertainty function 
ux) = X m (y); x, y € 2°’ (6) 
VOAX#D 


yx 
The sum of the belief, disbelief, and uncertainty 
values is equal 1: 
b(x) + d(x) tux) =1;xE€2% xð (T) 
As noted above, the set 2 can include both 
simple and compound elements-hypotheses. Let x 
be particular subset of set 2, then the number of 
elements it contains, is called atomicity of x and 
noted by ixl. With help of the atomicity concept, 
one can estimate the degree of overlap that exists 
between any two subsets of the set 2'!. 
Suppose x, y € 2, then for any given y # Ø the 
relative atomicity of x to y is defined as 


| 
le le xn yl 
y ly! 
Note that xay = Ø => a|) oan 
y 





See 2” yA. (8) 


x NPP 

xcy>a = 1. That’s, relative atomicity lies 
y 

within the [0, 1] interval. The relative atomicity of 

the atomic element (state) to its FoD, denoted by 


Xa : : 
a (=) , is written as a(x). Unless otherwise stated, 


the relative atomicity of an element (state) is 
referred to the FoD. 

The subjective logic also defines for given 
element x the probability expectation function 
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EQ) =>) me(y)-a =; xye2 9) 


Quite often, a researcher is not interested in all 
elements of the FoD @, but only in one particular 
element. In such cases, in order to simplify the 
calculations, it is proposed in subjective logic to 
derive from an initial FoD © a binary FoD 6” 
capable the following properties [15]: 

A FoD @* contains only two atomic elements: 
element of interest x and —x where —x is the 
complement of x in the FoD ©. 

The BMA m on binary FoD ale is defined 


according to: 
ms, (x) = D(x) 


mz (=x) = d(x) 
m, (0*) = u(x) 
A binary FoD Ø* with the BMA Mz. that 


satisfies (10) is called the focused FoD with focus 
on x. The focused relative atomicity of element x is 
defined by the following equation: 


E(x)—b(x 
= 
u(x) 
a;.(x)=a(x), if u(x) =0 
It should be emphasized that the belief, disbelief 
and uncertainty functions of x are identical in set 


(10) 


if u(x) #0 
(11) 


2'°' and set o’. The focused relative atomicity is 
also defined so that the probability expectation 
value for the element x is the same in the sets © and 
ð * 

The operation of transforming a FoD @ in a 
binary set Ø* is useful in situations when it is 
necessary to combine hypotheses. 


4 The method 

Our ultimate purpose is to develop the method 
for assessing the oil and gas potential of every plot 
within the given area. Structural and logical scheme 
of the proposed method is shown in Fig. 1. 
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Fig. 1. Structural and logical scheme of the proposed method 


The method is a multi-step one, the work begins 
with the collection of data on the study area and 
should end with the creating an information product 
in the form of a map with the division of the plots 
on classes in accordance with their forecasted oil 
and gas potential. 

Step 1. Data collection and preparation. The 
amount of data required is large enough, so usually 
all available technical and information sources are 
exploit, including: 

= equipment for seismic survey; 

=" magnetometric equipment; 

= electrical intelligence equipment; 

= aerial and satellite imagery; 

= geological survey data; 

= spectroradiometric devices; 

= humidity maps and data of soil geochemical 
analysis; 

= topographic maps and maps of lineaments; 

= digital elevation models, digital terrain models 
and many other. 

Using technical means, the necessary 
observations and measurements are carried out, and 
they are planned in such a way as to obtain 
sufficient information on the spatial distribution of 


E-ISSN: 2415-1521 


93 


the required physical parameters 
humidity, absorption, etc.). 

Some part of the required information is 
obtained through geological funds, by searching the 
global telecommunication networks, etc. The 
spatial distribution of physical parameter over 
study area is considered as a separate information 
layer and together all layers create a hypercube. 

It should be noted that the data obtained from 
different sources have different presentation 
structure, formats, differ in size, frequency of 
measurements, etc. [16]. Therefore, measures to 
ensure the interoperability of the collected data 
must be taken. 

Step 2. Data Analysis. Roughly speaking, all 
plots of study area can be divided on two 
categories: having hydrocarbon potential (so named 
positive) ones and ones don’t containing 
hydrocarbons (negative). As the positive plots are 
considered oil, natural gas and condensate 
collectors. 

Therefore, initially we have four hypotheses 
regarding possible state of a plot: 

A, — oil; 

Hh — natural gas; 


(brightness, 
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H3 — condensate; 
H; — there are no signs of hydrocarbons. 
These four hypotheses are simple ones and 


together create the FoD. 
The purpose of data analysis is to obtain for 
each plot hypotheses distribution and 


corresponding belief masses. Such distributions 
must be form separately and independently for each 
information layer. 

Suppose, there are K plots in the study area and 
as a result of the measurements and observations on 
step 1, and L information layers number have been 
constructed. Hypotheses for the k plot 
(k= 1, 2,...,K) obtaining by analysis of the data of 
the J layer (J = 1,2,..., L), will be indicated H'(k). 

The set of hypotheses regarding the state of k 
plot (k = 1, 2,...,K), which is formed on the basis of 
the analysis of the data of the / layer, can be written 
as a tuple 


Q, 
q=1 
where H : (k) is the q} hypothesis regarding the 


(H, (k);m; (k) (12) 


state of the k plot; m, (k) is the belief mass of the 


hypothesis H (k); Q, is a total number of 


hypotheses regarding the state of the k plot, which 
are formed on the basis of the data of the / 
information layer. 

It must be pointed out that the task of forming 
tuples (12) requires considerable intellectual effort; 
therefore, it is assumed that an expert method is 
used to solve it. In the set of hypotheses (12), 
which is formed for a plot, both simple and 
compound hypotheses can be. The next step is to 
form the opinion about the state of every plot of 
study area. 

Step 3. Conversion of hypotheses into opinions. 

Each simple hypothesis from the FoD is 
converted into the corresponding opinion. For this, 
based on the known belief masses of this 
hypothesis, the values of the belief, disbelief and 
uncertainty that determine this corresponding 
opinion must be calculated using expressions (4)- 
(6). 

Calculating the mathematical expectation of the 
hypothesis, it should be borne in mind that the 
correct transition from the initial set © to the binary 


set ©" is hindered by the interdependence between 
the probability expectation function E (x) and the 
focused relative atomicity a, ( x) _This problem 


may be solved with help the Hierarchical DSmP 
transformation procedure proposed in [17]. Its 
essence is as follows. 
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Suppose there is the FoD © = f H,,H,, H,} and a 


known the set 
2°’ = { (Hh; m),(H>; m2),(Hs; ms), 
(Hi2; m12),(133 m3), 
(H23; m23),(H123; M123) } 
in which each simple or compound hypothesis H 
has own corresponding belief mass m. Then the 
procedure of hierarchical probabilistic 


(13) 


transformation can be explained with help of Fig. 2 





Fig. 2. Graphic illustration of Hierarchical DSmP 
transformation 


Transformation process starts with the re- 
distribution of the belief mass of the compound 
hypothesis (Hı23; m123) = (Hı U H2 U H3; miz23} to 
the belief masses of the compound hypotheses 
(Hı U Hh; mı27, (Hı U H3; m3), (H2 U H3; m23). 
This stage occurs by adding to their masses the 
product of mass mı23 of the compound hypothesis 
Hi3 by the re-distribution coefficient, which 
characterizes the proportion of the corresponding 
mass in the total sum of masses H12, Hi3 and H23. 

In the second stage, the two-membered 
disjunctions are decomposed into singletons, and 
for each of them the belief mass value is 
determined in the same way. Re-distribution of 
masses between hypotheses in the given stages is 
carried out in accordance with the expression: 

ATpsn(A(v-1)) = m(H(v—1)) + 

+ )[m(H(v))- RHO -1)),A()) 


H(v)>H (v-1) 
H(v),H (v-l)e2™ 


where R(H(v-1),H(v)) is a proportional re- 
distribution coefficient, which is defined as 
m(H(v—1)) 
>, m(A(v-1)) 


H (v)>DH (v-1) 
H (v), H (v-1)e2® 


where v is a stage of transformation. 

The described procedure is implemented into an 
algorithm 1 represented below. The algorithm 
provides the formation of an opinion space; in 
pseudocode, it is assumed that the total number of 
simple hypotheses is J. The hypotheses are formed 
for every plot of the study area in each layer 


(14) 


R(A(v-1),A(v)) = (15) 
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separately and independently. Further, an opinion 


space is algorithmically formed from these 
hypotheses, 1.e., the algorithm performs the 
transformation 

How, =O d a a, N (16) 


E Doan ancl T a 


Algorithm 1: The algorithm to generate an opinion 


sets 
Input: Belief Mass Assignment 
// Form as 
m= ({ Him}, Em hal Asm, }, 
[=1,2,..,L,q=1,2,...,0 
for q in O': 

b(H;)= X, mo(H,) 
H,cH; 


d(H')= > m,(H!) 


l Be 
HOH j= 


ZS I 
u(H,)= 2 mo(H,) 
H OH ;4#0 
HiH! 


E(H\(v-1)) = m(H!(v-1)) + 
+ > [mE V) RH v -1), H v) 


H! (v)DH; (v-1) 


mH(v-1 
RH 0-1), H\(v)) = ae 
H (v)DH (v-1) 
H (v), H (v-l)e2 
if u( H`) =() 
a( H`) ="), 
else 
,,_ E(H;)-b(H;) 
a) 
j 
end 
Output: w = (w; , W}, a., w, pW 


l=1,2,...,L; j=1,2, ...,J 


It should be noted that the opinions w are 
determined by the belief masses of hypotheses H, 
which are estimated at step 2. Both estimates of the 
belief masses and decisions are usually carried out 
under conditions of partial uncertainty. Such 
conditions of data analysis adversely affect, in 
particular, the accuracy of the belief b and disbelief 
d values. 

It was shown in [15] that under such conditions, 
belief b and disbelief d functions acquire greater 
adequacy if the uncertainty fractions that are 
present in these functions are exported into the 
uncertainty function u. This can be realized by 
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adjusting the uncertainty function u values using 
the uncertainty maximization procedure. 

Step 4. Adjustment of belief and disbelief 
functions into opinions. The essence of the 
uncertainty maximization procedure is as follows 
[18]. Suppose that w; = (bj, dj, uj, aj) is a definite 
opinion about a binary unique event. Then at 
known b; and d; uncertainty maximization of an 
opinion w; is carried out by its transformation into 


an opinion w, = (b,,d.,u.,a,) as follows: 


for E(wj) < qj : for E(w;) > qj: 








= : d. 
ae a ae ree 
b. La, 
2.214 x 

a, d, =0 

5 (17) 

7 ee u =u, et —— 
J J a, J J —a, 
a, =a, a,=4a, 


Thus, for any k plot (k= 1, 2, ...,K), a set of L 
tuples of opinions (by the number of information 
layers) is formed: 








(wi (k)) i fa (Wi (ke). (ke)... (kD) 5 
(wi (k)) WP (Vif (k), Wi (k), (KD) 5 
(wi (k)) a (WK), (Kacey (K)) 18) 





The next task is to combine these L tuples into a 
single tuple whose components will reflect the 
properties of the k plot (k = 1, 2, ...,K). 

The algorithm for the uncertainty maximization 
procedure is presented below. 


Algorithm 2: The algorithm of uncertainty 
maximization process 
Input: opinion sets. 
// Form as w' = (w! , w}, «.. 5 W ness Wp), 
feo ye er! Bit a ea eN | 
for j in J: 
if E(w,) > a' 





Volume 8, 2020 


WSEAS TRANSACTIONS on COMPUTER RESEARCH 
DOI: 10.37394/232018.2020.8.11 





u.=u,+ 
J l 
la. 
Al] 
a,=a, 
else if E(w,) <a 
a 
b, =0 
b' 
d,=1-u,-— 
«i 
l 
a. =u" + — 
J a! 
j 
Al Oo l 
a,=a, 
end 
Output: w'=(w,,W,,..., W oy w; ); 
f=), 2; re Ok Eo eee 


Step 5. Combining opinions. Opinion tuples for 
each plot are combined independently. We assume 
that for any plot we have L tuples, each of which 
contains J opinions. Tuples are combined 
sequentially: first, the tuples of opinions formed 
from data of the first (‘=1) and second (l= 2) 
information layers are combined. Further, the new 
tuple, which is obtained as a result of combining, is 
combined with the tuple previously formed from 
the data of the third (l = 3) information layer, and 
so on. 

Consider the mathematical aspect of combining. 


Suppose there are two tuples: 
(Vii (K), (k) gees, (k)), 
( w” (k),... W” (k),... WEO), 


and into them the opinions of the same name must 
be combined in pairs ŵ' =(b', d', u,a and 
yw =(b"!, d™!, u l+1 va: 
In subjective logic, the operation of combining 
opinions is denoted by a symbol ® and is written as 


follows: 


w” (k) = Ww, (k)® w” (k) (19) 

The result of combining w" (k) = 
(b; d", u", a") is called a consensus 
opinion, and its belief, disbelief, uncertainty 


functions, and atomicity are calculated using the 
following formulas [18]: 
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l I+] 141_ 1 
bi = bu +b, u, 
J aa 1_ I 
u’ + u, —UuuUu. 
J J 
l I+] 1+1_ 1 
a = du +d, u, 
J [+1 1_ l+ 
u’ a Zuu, 
2 ae (20) 
LIH __ U Ul 
j E 7 1 l l+ 
" u, +u, “Hu i : 
l y” ES ae) l IHN. | 
mau 4U; ta; u, - (a, +a; )u, 
j E a l I+] 
i u’ +u + — 2u u : 


The pseudo-code in Algorithm 3 illustrates the 
fusion process intuitively. 


Algorithm 3: The fusion process 
// ® represents using consensus operator for 
combining opinions 
Input: opinion sets. 


// Form as 

We (Wi shi Wc Wy eas: 

WW, w o We jan We o 

VW, T e WD 
PT Dien pA Pe Mele eee ed 

for j in J 


W= ORO... DM OD... DN 


j 
end 


Output: W = (Wi, W2, ..., W;, ..., WJ); 
Pa E | 
Further, the obtained consensus opinion 


l 


(w (k)) d 


is combined with next opinion — 


the tuple of (18). In total the operation of 
combining opinions is performed L-1 times, and as 
a result one tuple containing up to J different 
opinions is calculated for the plot - 
W = (Wi, Wo, ..., Wj, ... , WD). 

In the same way, the tuples for all plots are 
calculated. 

As a rule, each tuple contains different 
consensus opinions about the oil and gas potential 
of the proper plot. It remains, to compare these 
opinions which each other and to choose one 
opinion that evaluates geological reality the most 
accurately. 

Step 6. Forming final opinions on gas and oil 
potential of subsoil plots of the study area. The 
final assessment for each plot is obtained by a 
sequential pairwise analysis of all calculated 
opinions on its oil and gas prospects (see tuples of 
opinions in (18)). The procedure for comparing a 
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pair of opinions and choosing the most preferred 
(most accurate) opinion is carried out in accordance 
with the following algorithm. 

Let there be a pair of opinions, one of which is 
W,, another is W2. The W; opinion is characterized 
by E(W;), u(Wi), a(Wı) values, W2 opinion is 
characterized by E(W2), u(W2), a(W2) values. 

The procedure for a comparative analysis of 
these opinions begins with checking the conditions: 

IE(W,) — E(W2)| < ££ (21) 
where sz = 0.1 min {E(W,), E(W2)}. 

If condition (21) is not satisfied, then choose is 
made in favor of the opinion characterized by a 
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bigger value of E(-). If condition (21) is satisfied, 
then next condition is checked: 

lu(W1) — u(W2)| < €n (22) 
where €, = 0.1 min {u(W,), u(W2) }. 

If condition (22) is not satisfied, then the 
opinion with the biggest uncertainty value u(-) is 
choosed. Otherwise, the values a(W,) and a( W2) are 
compared and the opinion with the biggest relative 
atomicity value a(-) is choosed. 

The described algorithm is sequentially applied 
to assess the oil and gas potential of each of the 
plots of the study area (Fig. 3). 
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Fig. 3. Flowchart of the algorithm for the most accurate opinions determining regarding oil and gas potential of 
plots of the study area 


Finally, according to the results of obtained 
unambiguous assessments of oil and gas potential 
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of the plots, a map of oil and gas content for all 
study area is constructed. 
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Loading the proposed method into a computer 
environment, one should take into account its quite 
big computational complexity. The greatest 
computational resources are required to implement 
Algorithm 3. The computational complexity of the 
matrix operations necessary for its implementation 
can be estimated using the expression O (n,m) ~ 
2*°xnxmxlogom, where m is the number of 
information layers, n is information layer 
parameter. Calculations show that for m=10; 
n=(255x257) ~ 2!° (these figures are taken from the 
practical example below) we have O (n,m) ~ 2* 
operations. 

In general, the proposed method belongs to the 
class of tasks having the exponential computation 
complexity; however, recently in [22] a highly 
efficient exact dynamic programming algorithm 
was proposed allowing successfully to solve similar 
tasks. 


5. Testing the proposed method 
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The developed method was tested under the 
Lypovodolyn oil and gas condensate field. In 
tectonic terms, the field is located in the western 
part of the paraxial zone of the Dnieper-Donetsk 
depression and is a part of the Artyukhiv- 
Lypovodolyn shaft. In 1962, the Lypovodolyn 
uplift was discovered by — structural-mapping 
drilling in Paleogene deposits. According to data of 
seismic works, structural maps of the covering of 
horizon IV, horizons Vb2 and Vc3 were 
constructed. In 1972, parametric well No. 453 was 
drilled in the area. Seismic studies identified two 
coves — Lypovodolyn and Penkiv [19]. 

Geomorphologically, the area of the 
Lypovodolyn oil and gas content structure (Fig. 4) 
is located on the right bank of the Khorol River 
within the Neo-Kharkiv Neogene terrace [20]. In 
Fig. 4 productive and nonproductive wells with 
numbers, major faults and lineaments, iso-gypsum 
of horizon B26 with depth markers, neotectonic 
uplifts, morphotectoisogypsum with elevation 
marks are shown. 
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Fig. 4. Lypovodolynsk oil and gas field 
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The coves of the structure 1s accentuated by the 
discordant laying of the right-hand tributary-moat 
of the river Horol, the top of which branches out, 
covering the raised relief block within the coves of 
the structure on both sides, forming a forked plan. 
This block is clearly fixed by the short left and 
right tributaries of the main moats flowing from the 
top of the block, which is emphasized by the plane 
flushing of the surface soils from the plane of the 
block and indicates its neotectonic activity. Similar 
in activity block is observed northwest of well 
No.9, from which it is separated by a lineament. 

The initial data for integration at the 
Lypovodolynsk deposit are the spatial layers 
contained in the corresponding files of various 
formats and coordinate systems: 

= Landsat/ETM+ satellite images, 6 spectral 

bands; 

= SRTM terrain v3 digital terrain elevation data; 

= land surface temperature; 

= lineament density raster maps, 4 pcs. in 45° 

increments; 

= vector map of structural horizons; 

= Neotectonic block map; 

= map of geochemical anomalies; 

= areas of optical anomalies of vegetation; 

= location of productive / non-productive wells; 

= map of depth of the terrain breakdown. 


The estimation of oil and gas potential by the 
developed method was carried out on the basis of 
data from 10 information layers (L=10). Tuples of 
primary data (see expression (12)) were formed by 
the expert method. Highly qualified geological 
experts were recruited; each of them had access to 
all information layers, was working independently 
and formed their own estimation hypotheses with 
assigning appropriate probabilities. 

The problem of estimating oil and gas potential 
of the area under study was solved by two methods: 
by a method based on the Dempster-Shafer theory 
of evidence and described in [12] and by the 
proposed method. 

The result of application of the method [12] for 
the forecasting oil and gas potential for different 
plots within the Lypovodolyn field is shown in Fig. 
5. From its analysis, it is clear that for a large area 
of the deposit the value of the belief function is 0. 

This is due, first of all, to the fact that, as shown 
in [21], applying Dempster’s combination rule in 
conflicting situations produces an unreasonable 
results, especially when belief masses of some 
hypotheses are close on zero. 

The result of application of the proposed 
method is shown in Fig. 6. From its analysis it is 
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clear that the entire area of the field is classified for 
the presence of hydrocarbons. Thus the proposed 
method is relieved of the disadvantages that are 
typical for Dempster’s combination rule and allows 
to obtain not a blurred interval but an unambiguous 
point estimate. 





= |o-01 E o2-03 |) 0.4-0.5 E o7-08 
E 01-02 E 03-04 D o5-07 M os-1 
Fig. 5. The result of using Dempster-Shafer method 
[12] for the Lypovodolyn field 


According to the results of the analysis of oil 
and gas content with the help of the proposed 
method the map of the distribution of output 
probabilities, by means of which it is possible to 
provide a complex evaluation of oil and gas 
potential of the studied area, has been plotted up 
(Fig. 6). 
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Fig. 6. The result of using the proposed method for 
the Lypovodolyn field 


It should be noted that using the proposed 
method emphasizes the block heterogeneity of the 
deposit better, which is confirmed by geological 
exploration. Unproductive well No. 10 is classified 
as high predictive probability because it is at the 
edge of the neotectonic uplift. As for productive 
wells, they all relate to 80-90% forecast probability, 
as evidenced by drilling data. 


6 Conclusions 


The new approach and method to forecasting 
the oil and gas potential of subsoil plots is 
proposed. This approach actively exploits the 
mathematical apparatus of subjective logic. 

The advantages of the proposed method are that 
it allows to adjusting the subjectivity and 
uncertainty of expert judgments, to generate 
adequate decisions on the results of the analysis 
one-time or rare events, and to avoid the majority 
of logical-computational problems that usually 
arise when combining conflicting hypothesis and 
opinions. 

The proposed method was tested under the 
Lypovodolyn oil and gas condensate field. 
According to the results, the proposed method 
showed advantages over the method based on the 
Dempster-Shafer theory of evidence. 

It is assumed that further research and efforts of 
authors will be focused on the technologization of 
the proposed method, the improvement of its 
algorithmic support and the writing of 
recommendations for its using in oil and gas 
exploration works. 
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